A critical analysis of vacancy-induced magnetism in mono and bilayer graphene 
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The observation of intrinsic magnetic order in graphene and graphene-based materials relies on the 
formation of magnetic moments and a sufficiently strong mutual interaction. Vacancies are arguably 
considered the primary source of magnetic moments. Here we present an in-depth density functional 
theory study of the spin-resolved electronic structure of (monoatomic) vacancies in graphene and 
bilayer graphene. We use two different methodologies: supercell calculations with the SIESTA 
code and cluster-embedded calculations with the ALACANT package. Our results are conclusive: 
The vacancy-induced extended n magnetic moments, which present long-range interactions and 
are capable of magnetic ordering, vanish at any experimentally relevant vacancy concentration. 
This holds for cr-bond passivated and un-passivated reconstructed vacancies, although, for the un- 
passivated ones, the disappearance of the n magnetic moments is accompanied by a very large 
magnetic susceptibility. Only for the unlikely case of a full a-bond passivation, preventing the 
reconstruction of the vacancy, a full value of \\ib for the 7r extended magnetic moment is recovered 
for both mono and bilayer cases. Our results put on hold claims of vacancy-induced ferromagnetic 
or antiferromagnetic order in graphene-based systems, while still leaving the door open to cr-type 
paramagnet ism . 

PACS numbers: 72.80.Vp,73.22.Pr,72.15.Rn 



I. INTRODUCTION 

According to theory, the existence of intrinsic (with- 
out invoking foreign species) magnetism in graphene and 
graphene-based materials should be the rule rather than 
the exception. Aside trivial paramagnetism associated 
with a dangling bonds of undercoordinatcd C atoms, 
a more interesting 7r magnetism should arise at zigzag 
edgesir— , bulk defects such as vacancies^ - — , or grain 
boundaries 2 ^, which either appear naturally or can be 
created in a more or less controlled manner. However, 
undisputed experimental evidence of magnetic order re- 
mains elusive since early claims of observation of ferro- 
magnetism in irradiated graphit o 21 ' 22 and graphene 2 ^. 
Recent claims based on transport^ and local probe 
measurement o 20 i 25 ' 26 seem more solid, but not entirely 
free from controversy 2 ^. The reasons why the observation 
of ferromagnetism in graphene and graphene derivatives 
is so elusive, even at low temperatures, are still unclear, 
but can generically be traced back to two facts: 1) the 
magnetic instability leading to the appearance of mag- 
netic moments can be superseded by structural (Jahn- 
Teller) instabilities or unwanted passivation by foreign 
species, and 2) the underlying antiferromagnetic corre- 
lations inherent to graphene favor this type of magnetic 
order over ferromagnetism even if the magnetic moments 
truly exist. 

Graphene represents the paradigm of bipartite lattices. 
At the heart of the bipartite nature lies the reason why 
some graphene derivatives result in half-filled 7r states at 
the Fermi energy which spin-split due to electron-electron 
interactions. When these interactions are restricted to 
be local, as described, e.g., by a one-orbital Hubbard 
model, and the electron-hole symmetry is exactly pre- 
served, the existence of a magnetic ground state with to- 




FIG. 1. (Color online) Atomic structure around a single atom 
vacancy in graphene. The arrow indicates the atom with the 
dangling a bond and the colored atoms are those rebonded. 



tal spin S — \Na — Nb\/2 is guaranteed by a theorem by 
Lieb 28 , where Na—Nb is the difference between the num- 
ber of atoms in each sublattice, i.e., the sublattice imbal- 
ance. This imbalance appears whenever graphene is cut 
or grown into triangular shapes bounded by zigzag edges 7 
or, conversely, when C atoms are removed from bulk 
graphene, creating vacancies^ or voidsi^ with similar tri- 
angular shapes. Interestingly, even when Na — Nb = 0, 
as is the case in zigzag nanoribbonai~— , large hexago- 
nal graphene nanoflakes with zigzag edgea 7 -^, or voids 
of similar shape in bulk graphene^, magnetic solutions 
may appear, but always with an envelope antiferromag- 
netic order on top of a local ferromagnetic order. Smaller 
structures with Na — Nb = are not magnetic as, e.g., 
recent work on divacancies in graphene 2 ^ shows. 

In between Lieb's theorem and the observation of mag- 
netic fingerprints in graphene-based systems stands a 
number of assumptions easily overlooked, namely; i) that 
the bipartite atomic structure is preserved, ii) that hy- 



drogenation or passivation of the extended ir states does 
not occur, iii) that the Hubbard model is a good ap- 
proximation to describe graphene 7r interactions, and 
(iv) that the substrate does not play a significant role. 
Despite that both chemical synthesis and physical ap- 
proaches are making progress into the creation of locally 
or globally imbalanced graphene structures, it still re- 
mains a remarkably challenging task. The second condi- 
tion is a major experimental challenge since it is difficult 
to avoid full passivation of the edges^ (not to mention 
to achieve the often implicit selective passivation of the a 
bonds, while avoiding that of the ir states), at least under 
conditions compatible with standard magnetic measure- 
ments. Moreover, if passivation is completely avoided 
in ultra-high vacuum conditions, the equilibrium atomic 
structure may develop a reconstruction that ruins the 
first condition and that, energy wise, competes favorably 
with the magnetic instabilit y 31 ' 32 . Substrates may pre- 
vent this reconstruction from taking placed 3 -, but it is un- 
clear whether or not they always respect the magnetic in- 
stability. Finally, while the (mean-field) Hubbard model 
has shown its reliability in reproducing results obtained 
with more sophisticated approximations [typically den- 
sity functional theory (DFT)£] , the comparison has only 
been carried out in the most favorable situation, namely, 
saturating the a bonds with H, thus avoiding unwanted 
lattice reconstructions. The extent to which the unsat- 
urated a bonds or lattice reconstructions may invalidate 
the use of the one-orbital Hubbard model remains largely 
unexplored. 

While the magnetic instability at zigazg edges is under 
present theoretical and experimental scrutinity, the nat- 
ural hosts of magnetism in bulk, vacancies, have not re- 
ceived due critical attention so far. Single-atom vacancies 
in bulk graphene are the simplest structures complying, 
in principle, with the conditions for the appearance of 
extended ir magnetic moments, both in mono^ ' 16 ' 18 ' 19 ' 34 
and multilayer graphene 3 ^—. When H saturation of the 
dangling bonds (left upon removal of a C atom) prevents 
the Jahn- Teller reconstruction of the vacancy, the value 
of the induced magnetic moment is expected to be 1 /is 
for the it electrons plus or 1 fis for the a bonds, depend- 
ing on whether 3 or 2 H atoms are available for satura- 
tion, respectively. Discrepancies, however, can be found 
in the literature regarding the actual value of the mag- 
netic moment when a single H (or no H at all) saturates a 
dangling bond and the vacancy reconstructs (see Fig. [1}. 
Values for the tt magnetic moment ranging from f=a 0.0 [is 
to « 1.0 fiB have been reported in this cas o 11 ' 12 ' 14 ' 16 ' 19 ' 38 . 

To illustrate the source of the discrepancy we show in 
Fig. [2] the electronic structure of a single vacancy in a 
graphene monolayer in two different instances. In panel 
(a) we consider the three a dangling bonds, left by the 
removal of a C atom, saturated with H atoms, whereas 
in panel (b) we consider no H passivation. In both cases 
a full atomic relaxation is carried out. A schematics of 
the atomic structure in the second case is shown in Fig. 
[T] which is similar to the ones previously reported in the 
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FIG. 2. Calculated spin resolved band structure (left panels) 
and densities of states (right panels) at a vacancy in a 6x6 
super-cell in monolayer graphene. In (a) the three a dangling 
bonds left by the removal of C atom are saturated with H 
atoms. In (b) a bare vacancy is considered. Full relaxation 
of the atomic geometry is performed in both cases. Solid 
and broken lines indicate spin majority and spin minority 
electronic states, respectively. The zero of energy is at the 
Fermi level. 



literatur o 11 ' 12 ' 14 ' 16 ' 19 ' 38 . The differences in the electronic 
structure between both situations are remarkable. In the 
former case the trigonal atomic symmetry is maintained 
and, despite small deviations from the perfect atomic lat- 
tice, Lieb's theorem is still expected to apply. We obtain 
two well separated spin minority and spin majority 7r 
bands (dispersive due to the supercell periodicity) near 
the Fermi level situated at the Dirac point. The mag- 
netic moment associated with the 7r orbitals is actually 
1 hb- In Fig. Hfb), aside the appearance of a er band 
at -0.75 eV, the situation is different. The it bands over- 
lap in energy close to the Fermi level and, therefore, the 
magnetic moment is smaller than 2 /is (1-71 in this case). 
Importantly, the Dirac point lies above the Fermi level 
at around 0.25 eV. 

These results already show the subtle effect of breaking 
the bipartite character of the lattice by the relaxation of 
the atoms around the vacancy which, in turns, induces 



the breaking of the electron-hole symmetry. In the case 
of bilayer graphene the situation is expected to be even 
more subtle since, along with these details, the influence 
of the bottom layer needs to be taken into consideration. 
Graphene layers are coupled via van der Waals interac- 
tions and, to our knowledge, the existing calculations did 
not tackle this issue appropriately. Also, as discussed be- 
low, the value of the ir magnetic moment turns out to be 
quite sensitive to details of the calculation such as size 
and shape of the supercell, generally decreasing with size. 
To compound things even further, a recent scanning tun- 
neling spectroscopy (STS) observation of the density of 
states associated with vacancies created under controlled 
conditions on the surface layer of graphite has been in- 
terpreted as a manifestation of ir magnetism 25 . A second 
reading of the experimental data, however, can also be 
interpreted otherwise, as the evidence of the absence of 
magnetism. This has prompted us to question whether 
vacancy-induced 7r magnetism exists at all. 

Here we present extensive and detailed DFT calcula- 
tions for vacancies in both mono and bilayer graphene, 
including van der Waals interaction^ as implemented 
in the SIESTA codei^i 2 -. Our results indicate that re- 
constructed vacancies in mono and bilayer graphene can 
host highly localized a magnetic moments of 1 /is (i.e., 
one unpaired spin) , but that the overall extended n mag- 
netism progressively decreases as the size of the supercell 
increases, extrapolating to zero in the zero- concentration 
limit. This limit is reached more rapidly when the dan- 
gling a bond is saturated with H (and the a magnetic 
moment is quenched) and also when the vacancy is cre- 
ated on the bilayer. In the latter no significant differ- 
ences are appreciated between the two lattice sites in 
this regard. Calculations for a single vacancy performed 
with a cluster-embedded methodology as implemented in 
the ANT.G code^ also yield values of the magnetic mo- 
ment in the 7r orbitals approaching zero, strengthening 
our conclusion. On the other hand, as expected, a value 
of 1 /is for the it magnetic moment is obtained for the 
unlikely case of a total hydrogenation of the a orbitals 
which prevents the Jahn- Teller reconstruction and recov- 
ers the Lieb's scenario. 



II. METHODOLOGY 

Most of the calculations reported here have been per- 
formed with the SIESTA code which uses a basis of nu- 
merical atomic orbitals^ and separable^ 5 - norm conserv- 
ing pseudopotentials 4 ^ with partial core corrections^. 
We have found satisfactory the standard double-^ basis 
with polarization orbitals (DZP) which has been used 
throughout this work. Also a ghost atom at the va- 
cancy has been included to improve the basis set, al- 
though it does not change the DZP results. The conver- 
gence of the relevant precision parameters was carefully 
checked. The real space integration grid had a cut-off 
of 500 Ryd. Of the order of up to 600 k points were 



used in the two-dimensional Brillouin zone sampling us- 
ing the Monkhorst-Pack k-points sampling. Spin resolved 
calculations are performed in most cases. To accelerate 
the self-consistency convergence, a polynomial broaden- 
ing of the energy levels was performed using the method 
of Methfessel and Paxton^ which is very suitable for sys- 
tems with a large variation of the density of states in the 
vicinity of the Fermi level as is the case in our system 
(see below). Broadening like Fermi-Dirac can be inap- 
propriate and give wrong results. It is worth mentioning 
that the energy differences between non-magnetic and 
magnetic solutions are, in general, small, what requires 
a very high convergence in all precision parameters and 
tolerances. To obtain the equilibrium geometry we re- 
laxed all the atoms until the forces acting on them were 
smaller than 0.01 eV/A. We obtain for the defect free 
graphene layer a nearest-neighbor distance of 1.435 A as 
compared with the experimental value of 1.42 A. In the 
bilayer calculation including van der Waals forces the dis- 
tance between planes is 3.42 A whereas the experimental 
one for graphite is 3.35 A. To calculate the geometrical 
and electronic structure of defects, we use the supercell 
calculation method with nxm cells containing the defect 
for n and m integers and standard unit cell vectors. As 
a general comment concerning the geometry of the va- 
cancy, we obtain results similar to those reported in the 
literature; the structure remains planar, two dangling a 
orbitals rebond in a new weak bond of whose length de- 
pends on the super cell size, ranging from 2.05 A in the 
6x6 case to 1.93 A in the 15 x 15 one. The third a orbital 
remains non-bonded. In Fig. [T]we show a reconstructed 
vacancy in the middle of a 5x5 supercell which, in the 
calculations, is periodically repeated in two dimensions. 



III. MONOLAYER GRAPHENE 

Although we are interested in the case of an isolated 
vacancy, we are forced to consider a finite concentration 
of vacancies on the same sublattice; this is what one actu- 
ally does in supercell calculations when using electronic 
structure codes such as SIESTA. As already briefly dis- 
cussed in the introduction, we have first carried out a 
6X6 supercell calculation in two cases: i) with H atoms 
saturating the a dangling bonds [Fig. HJa)] and ii) with 
no extra H atoms [Fig. [2Jb)]. In the former case, which 
is not so relevant from an experimental point of view, the 
full passivation of the three dangling bonds almost com- 
pletely prevents the reconstruction of the lattice while 
in the latter a strong Jahn- Teller distortion takes place 
(see [1]) . Left panels show the band structure and right 
panels the total density of states (DOS). Only 7r bands 
are visible in Fig. Uta) whereas the unsaturated a dan- 
gling bond forms a band at « -0.75 eV in Fig. HJb). 
This band presents a large spin splitting and is almost 
flat as corresponds to a highly localized state. On the 
other hand the splitting of the n band is much smaller in 
both cases and presents a visible dispersion which is due 
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FIG. 3. Calculated spin resolved band structure at a vacancy 
in increasingly nxn supercells in monolayer graphene. Panels 
(a), (b) and (c) stand for 9x9, 12 x 12 and 15 x 15 supercells 
respectively. Solid and broken lines indicate spin majority and 
spin minority electronic states. The zero of energy is at the 
Fermi level. The arrows in panel (a) indicate states induced 
by the defect [similar bands can be identified in panels (b) 
and (c)]. 



to the always present interaction between vacancies due 
to the semi-localized character of the ■n state created by 
the vacancy^. In the DOS the majority spin a peak and 
the two spin-resolved peaks coming from the 7r state are 
visible. 

For the fully saturated vacancy neither the unoccupied 
band nor the occupied one crosses the Fermi level. The 
magnetic moment is thus always quantized to 1 \lb (we 
have checked that this is case for any concentration of va- 
cancies), as predicted by Lieb's theorem. Note that the a 
bonds are saturated and do not host any magnetism here. 
Note also that the bipartite nature of the lattice is pre- 
served, essentially restoring the electron- hole symmetry. 
The low-energy physics resulting from these type of va- 




1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 
Magnetic Moment (u B ) 

FIG. 4. (Color online) Total energy versus magnetic moment 
for different supercell sizes for a monovacancy on a graphene 
monolayer. The inset indicates the magnetic moment at the 
total energy minimum as a function of the vacancy concen- 
tration (inverse of the supercell size). 



cancies is completely equivalent to the physics of hydro- 
genated graphene^ and has been discussed at length in 
Ref. [l5|. For the unsaturated vacancy the band structure 
presents subtle differences. Both spin-split 7r bands cross 
the Fermi level, the upper one actually staying pinned 
to it. This prevents the magnetic moment from reaching 
the saturation value of 2/ib (1 fJ-B from the a bond plus 
1 fj,B from the vacancy- induced it state), yielding a total 
value of around 1.71 [is for this particular calculation. 
This is linked to the remarkable fact that the Fermi level 
lies below the Dirac point, which is equivalent to saying 
that the vacancy acts as an acceptor impurity. 

We should note at this point that the value of the 
magnetic moment for the reconstructed vacancy, which 
is the relevant case from the experimental point of view, 
changes with the size of the supercell so we set out now 
to do a systematic study. Figure [3] shows the band struc- 
ture for an increasing supercell size sequence 3n x 3n up 
to 15 x 15. While the a band becomes quickly completely 
flat at around -0.8 eV, the 7r band retains the dispersion 
and the spin splitting although these become flatter and 
smaller, respectively, as the supercell size increases. This 
reflects the increasing distance between vacancies and the 
concomitant increasing extension of the it state induced 
by the vacancy at the Dirac point. (The lattice recon- 
struction does not allow us to establish a perfect analogy 
with the Dirac state in the standard tight-binding model 
which decays as 1/r, but we have no reason to expect 
otherwise). Interestingly, the partially occupied upper 
spin-split 7T band stays pinned at the Fermi level on a 
part of the Brillouin zone for all supercell sizes. Also 
the difference between the Dirac point and the Fermi 
level decreases which can be easily understood since the 
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FIG. 5. Band structure when the a dangling orbital of the 
vacancy is saturated with a H atom for a (a) 6 x 6 and a (b) 
12x12 supercell. Solid and broken lines indicate spin majority 
and spin minority electronic states. The zero of energy is at 
the Fermi level. 



concentration of vacancies (or acceptor impurities) de- 
creases. It is important at this point to notice that the 
vacancy does not only induce one 7r band around the 
Fermi energy. As indicated by arrows in Fig. [3] (a), a 
new set of 7r bands emerge which carry spectral weight 
of the resonance mostly in the occupied part of the spec- 
trum. In the zero-concentration limit, this set of bands 
should merge into a continuum and give a finite a width 
to the resonance in the energy sector of occupied states 
(see below). 

In Fig. @]we show total energy calculations as a func- 
tion of fixed magnetic moment fi for various supercell 
sizes. One can easily appreciate how, as the supercell 
size increases, the minimum energy value of /i, [1q, moves 
towards 1/^s, which is the lower limit imposed by the 
unpaired electron of the a dangling bond. At the same 
time, d 2 E(n)/diJL 2 \ tl0 — > 0, which amounts to a very large 



susceptibility per vacancy. This shallow variation of the 
energy with the magnetic moment is a remarkable fact; it 
should be noticed that, for instance, in the 15 x 15 case, 
the magnetic moment in the it states can vary around 
0.4 /j,b within 1 meV. This indicates that, even at low 
temperatures, the magnetic moment is ill defined. This 
is even more pronounced in the bilayer case (see below). 
The dependence of /iq on the inverse of the supercell size 
(i.e. the concentration of impurities) is plotted in the 
inset of Fig. 2] (see below for further analysis) . 

We finally examine the possibility of having the dan- 
gling a bond saturated with atomic H. The calculations 
are performed allowing relaxations of all the atoms as 
indicated above. Now the a band disappears from the 
energy window of interest (see Fig. [5]) along with the as- 
sociated magnetic moment. While for small supercells (or 
high concentrations of vacancies) the spin splitting of the 
7r band is still visible, it already completely vanishes for 
supercell sizes as those considered in the previous case. 
This result reflects the importance of considering the mu- 
tual influence between the a and 7r electrons, at least as 
far as magnetic properties is concerned, when the former 
are not part of a bond to other species such as, e.g., H. 
This effect cannot be captured by the Hubbard model 
where the saturation of the sigma bonds is always im- 
plied even if the hopping terms are adapted to the atomic 
reconstruction 51 . 



IV. BILAYER GRAPHENE 

We now consider vacancies on bilayer graphene with 
Bernal stacking. In this situation, removing a C atom 
from one sublattice or the other is different due to the 
underlying graphene layer, resulting in two types of va- 
cancies, a and f3, depending on whether or not the va- 
cancy is created on top or hollow position^ 5 -. The band 
structure of a vacancy in a 9 x 9 supercell in both cases 
is shown in Fig. [5] As in the previous case we allow 
for full relaxation of the atomic coordinates on the layer 
containing the vacancy. Figure [B](c) shows the bands for 
a bilayer without vacancies. The mass acquired by the 
Dirac electrons (the parabolic dispersion at the Fermi en- 
ergy) as a result of the interlayer interactions is evident 
in the plot. The tt bands associated with the vacancy, re- 
gardless of the sublattice creation site, are spin-split for 
small cells. Contrary to the monolayer case, the minor- 
ity spin band crosses the Fermi level, already indicating 
a stronger tendency towards the quenching of it mag- 
netism than in the monolayer. Therefore, it should not 
come as a surprise that, as in the monolayer case, the 
spin splitting goes to zero as the distance between va- 
cancies increases, remaining only the magnetic moment 
associated with the a bond. The difference between the 
a and [3 cases is minor. The bands in the a case are nar- 
rower than those of the (3 case as expected^ 5 .. The inset 
in Fig. [7] shows (1q as a function of the inverse of the su- 
percell size for a vacancies. Compared to the monolayer 
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FIG. 6. Band structure for vacancies on a bilayer for a 9 x 9 
supercell. Solid and broken lines indicate spin majority and 
spin minority electronic states. The zero of energy is at the 
Fermi level. Panels (a) and (b) refer to a vacancy created 
on top (a) or hollow (f3) positions, respectively. Panel (c) 
represents the defect free graphene bilayer. 



result in Fig. [H one can safely extrapolate /in — > 1/is in 
the zero-concentration limit. Note that the somewhat er- 
ratic behavior of /in as a function of the inverse supercell 
size can be attributed to considering all consecutive sizes 
while in the monolayer case we are only plotting results 
for 3nx3n supercells. In the light of the results, as for the 
monolayer, we also expect /xn — > 0/jb if the a dangling 
bond is saturated. Also, as in the monolayer case, we 
obtain shallow energy curves versus magnetic moment, 
indicating an even higher susceptibility (see Fig. [7]). 
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FIG. 7. (Color online) Total energy versus magnetic moment 
for different supercell sizes for an a vacancy in a graphene 
bilayer. The inset shows the magnetic moment at the total 
energy minimum as a function of the vacancy concentration 
(inverse of the supercell size). 



CRITICAL ANALYSIS OF RESULTS AND 
FINAL CONSIDERATIONS 



While the results for the bilayer vacancy and the mono- 
layer vacancy with H seem conclusive regarding the van- 
ishing value of the ir magnetism in the low concentration 
limit, the ones for the H free monolayer vacancy remain 
less clear. We would like to discard any possible influ- 
ence on the results of the specific sequence of supercells 
considered in our calculations so we have also performed 
additional calculations with supercells out of the main 
sequence in x 3n. We now plot in Fig. [8] all the results 
for the total magnetic moment, including values obtained 
with all types of supercells. In the light of this plot we 
can safely conclude that /j,n goes to 1 ub in the zero- 
concentration limit n — > 0, posibly as oc n where S < 1. 
In fact, despite that the behavior of u Q is not mono- 
tonic, a good fit to /in (n) = 1 + an + by/n can be done 
(o= -15.11,6=7.64). 

We would like to address now the influence of the 
periodicity on the results. To this aim, we have also 
performed calculations for truly isolated vacancies with 
the help of the ALACANT package^; in particular, we 
have employed our code ANT.G which interfaces with 
GAUSSIAN09 52 . In this case the supercell is sorrounded 
by an effective medium defined by a two-dimensional 
Bethe lattice^ of coordination three and Slater-Koster 
parameters for the C sp orbitals . Here the Green's func- 
tion of the supercell is selfconsistently computed subject 
to a fixed selfenergy representing the Bethe lattice. In 
contrast to the calculations with SIESTA, the vacancy is 
here trully isolated, but the electronic structure outside 
the cell remains fixed and unmagnetized. Unlike bulk 
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FIG. 8. (Color on line). Calculated magnetic moment induced 
by a vacancy in a graphene monolayer for various concentra- 
tions (inverse of the supercell size). The red circles are the 
calculated values and the broken line is a fit to 1 + an + b\fn. 



graphene, the Bethe lattice model presents a finite den- 
sity of states at the Fermi energy, which gives the quasi- 
localized 7r state of the vacancy a finite lifetime for any 
cell size even at zero energy. We have also used here the 
generalized gradient approximation through the BPBE 
functional as implemented in GAUSSIAN09^ and a ba- 
sis set equivalent to that in the SIESTA calculations. The 
atomic structure has also been optimized, obtaining es- 
sentially the same geometry. The values of the magnetic 
moments so obtained are all in the range ~ 1.1 — 1.3/is, 
with a clear trend towards 1.0/ib as the system size in- 
creases. One may conclude that the periodicity, if any- 
thing, enhances the values of the it magnetic moments 
induced by the vacancy. 

To make connection with available experimental 
information^, we have plotted in Fig. |H] the DOS pro- 
jected on the it orbital of the vacancy atom with the 
dangling bond (for the other two vacancy atoms the re- 
sults are similar). The calculation is a non spin resolved 
one for a 18 x 18 supercell. We obtain an asymmetric 
and almost fully occupied sharp-peaked resonance at the 
Fermi level, its spectral shape strongly deviating from 
a symmetric Breit-Wigner or 1/ \E\ resonance^. Most 
of its weight is in the valence band with no extra struc- 
ture in the conduction band and a small gap right above 
the main peak. This anomalous form of the line shape 
is a dramatic consequence of the electron-hole symme- 
try breaking (see qualitatively similar results in a model 
calculation by Pereira et al42). The asymmetry and the 
presence of the small gap right above the sharp peak 
would prevent, in the isolated vacancy limit, the Stoner 
instability and the formation of an extended magnetic 
moment. 

From these results several conclusions can be extracted 
regarding various experimental observations: 




FIG. 9. Density of states projected on the n orbital at an 
atom in the vacancy (solid line) in the non-magnetic solution. 
A small (0.05 eV) gaussian broadening has been included for 
presentation purposes. The corresponding density of states in 
a defect free graphene is represented by the broken line. 



Our results indicate that only a high concentration 
of ordered vacancies on the same sublattice can sus- 
tain finite values of the ir magnetic moments and 
lead to a ferromagnetically ordered state. The con- 
centration below which these magnetic moments 
disappear depends on whether or not the a dan- 
gling bond is passivated, being much higher for the 
passivated case. In addition, one should not for- 
get that, in average, the same number of vacancies 
are expected on both sublattices. In this case the 
7T magnetic moments are quenched when vacancies 
are in proximity^, further disfavoring the existence 
of these magnetic moments. On top of that, an ex- 
cessive concentration of vacancies will likely render 
graphene unstable. 

Although one should keep in mind that the STS re- 
sults by Brihuega et al<2£ refer to surface graphite, 
our results are compatible with their observations 
without invoking the existence of 7r magnetism. In 
their experiment no trace of two spin-split peaks 
near the Fermi energy can be seen. Furthermore, 
although the DOS in Fig. [9] corresponds to a 
graphene monolayer, the asymmetry in the exper- 
imental dl / dV peak at low bias nicely compares 
with our result. We should note, nevertheless, that 
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we obtain a large magnetic susceptibility mainly 
associated to the soft position of the spin-majority 
peak in the DOS. The possibility for thermal fluctu- 
ations to wash out this peak from the DOS, mask- 
ing the spin-split structure cannot be entirely ruled 
out. 

• As shown in Fig. [3] the a band becomes rapidly flat 
as the concentration of vacancies decreases. This 
indicates that these localized a magnetic moments 
do not interact for any reasonable concentration 
and should behave as paramagnetic centers. Upon 
completion of this work an experiment by Nair et 
al^ has unambiguously shown the paramagnetic 
behavior of irradiated graphene, ruling out any pos- 
sibility of magnetic order induced by vacancies and 
in complete agreement with our results. 

• To conclude, one should keep in mind that it mag- 



netism and magnetic order can still emerge through 
atomic H adsorption or through any other adsor- 
bate capable of similar covalent bonding to p z or- 
bitals. This magnetism should be amenable to ex- 
perimental verification, for instance in magneto- 
transport measurements, as recently propose d 55 ' 56 . 
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